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Abstract. Derivation of an analytic, closed-form solution for Q-switched laser induced op- 
toacoustic phenomenon in absorbing fluid media is presented. The solution assumes spherical 
symmetry as well for the forcing function, which represents heat deposition from Q-switched 
lasers. The Green's solution provided is a suitable kernel to generate more complex solutions 
arising in optoacoustics, optoacoustic spectroscopy, photoacoustic and photothermal problems. 
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1 INTRODUCTION 

Optoacoustic phenomenon is the formation of sound waves in absorbing media by means of 
either a modulated or pulsed light source OJ0. Regardless of the impulse, an acoustic wave 
involves the physical vibration of the constituent parts comprising the material and the cumu- 
lative propagation of the wave enjoins all material properties in total; as a consequence, the 
information contained within the wave envelopes all the fundamental physical properties of the 
material. The difficulty rest in deciphering signals detected by various means, where some 
detection methods measure the physical vibration and other methods exploit derivative thermo- 
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dynamic relations. 
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Photoacoustic and photothermal spectroscopy are two separate detection methods for the 
pi ' same physical phenomenon, that is, both are semantic nomenclatures for optoacoustic phe- 

nomenon. The generation of an acoustic pressure wave in ponderable media can be detected 
either by measuring the amplitude of the stress vector (the physical displacement) or indirectly 
by measuring the degree by which some intrinsic physical parameter changes. Direct measure- 
ment of the stress vector can be monitored by utilizing an electric transducer and this detection 
method constitutes photoacoustic phenomenon. Although many intrinsic characteristics change 
during the disturbance, measuring the degree by which the refractive index changes constitutes 
the photothermal detection method. This method makes use of the refractive index change 
caused by the pressure wave to deflect a probe laser beam. The probe beam is projected onto a 
differential photodetector, which differences the magnitude of several segmented regions of the 
photodetector. The probe beam deflects as the acoustic wave passes at point in space; further- 
more, both the direction and magnitude of the deflection angle is proportional to the shape of 
the refractive index change caused by the passing pressure wave. 

Laser technology enables ultra-short pulses of laser light, of which Q-switched lasers are 
an example. The physical interaction of photon energy and ponderable media is quite complex. 
The particular mode of interaction is dependent upon the wavelength of light, pulse energy, 
beam characteristics, pulse duration and material constants, such as thermal conductivity, den- 
sity, &c. Despite the many complexities of interaction, Lambert's absorption coefficient is an 
empirical measure of the coupling between light of a given frequency and media of interest. For 
high frequency radiation, the general mode of interaction is to excite both atomic and molecular 
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vibrational and/or electronic transitions. Excess energy absorbed by atoms and molecules are 
expelled through molecular relaxation processes. 

The optoacoustic phenomenon is an isentropic process, where the condensation (compres- 
sion) of the media at the leading edge of the wavefront causes the material to rapidly increase 
in temperature; similarly, as the pressure wave passes, the subsequent rarefaction (expansion) 
allows for the cooling of the material. This would be a reversible adiabatic process if it were 
not for nonconservative forces dissipating the energy as the acoustic wave propagates; albeit, 
the convolution of such nonconservative forces temporally and spatially allow atomic and/or 
molecular characteristics to be imprinted onto the resultant acoustic wave, which are then de- 
tected, identified and deciphered. Deconvolution of the resulting signal is confounded by the 
propagation of the disturbance through the volume of interest. Subsequent reflections from the 
container walls, material interaction and nonconservative forces from chemical constituents in- 
creases the complexity of the resultant wavefront; thus, detection of the initial disturbance is 
most efficacious. 

Optoacoustics and optoacoustic spectroscopy is undoubtedly a very rich and complex tech- 
nique and cannot be discussed comprehensively in such few words; but, the purpose of the 
present expose is to derive a compact, closed-form analytic solution for the coupled partial 
differential system describing this complex phenomenon. By compact is meant an analytic ex- 
pression relating all requisite parameters involved in generating an acoustic wave from a pulsed 
laser source in as few terms as possible. Thus, the intent is to present a functional that can serve 
as a kernel solution to generate more complex solutions, to be used in circumstances involving 
greater specificity in geometry, material properties, &c. Discussions will concentrate on the 
generation and behavior of the resulting acoustic wave; the topic of detection is deferred. 

2 ASSEMBLY 

For incompressible irrotational fluid flow, pressure p is proportional to the time rate of change 
of a scalar potential $ which describes the physical displacement of particles and fluid, encap- 
sulated by equation (1) (HE)- The scalar potential is related to the time rate of change of the 
temperature T in the fluid through the wave equation, where the characteristic speed of a result- 
ing acoustic disturbance through the medium is represented by v, enjoined by equation (2). The 
origin of any temperature disturbance results from a forcing function (S) acting on the fluid, 
defined by equation (3). Tracing the relationship the scalar potential has with pressure, it can 
be seen that pressure is proportional to the second-order time derivative of temperature within 
the fluid, hence, pressure is directly related to the particle and/or fluid acceleration. Tempera- 
ture is proportional to the absorbed portion of the laser's emitted time averaged radiant energy 
per volume; equivalently, the temperature is proportional to the absorbed portion of the emitted 
spectral irradiance [Watts per volume]. The temperature in the fluid is represented by the stan- 
dard forced heat equation, i.e. equation (3). The coupled partial differential system of equations, 
equations (1-3), model pressure waves initiated by thermal deposition from a forcing function 
acting on the fluid. 
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P° C p^T =«V 2 T + S (3) 
dt 

The forced heat equation describes the time rate of change in temperature, scaled by the fluid 
density (po) and the heat capacity (C p ), to be equal to the sum of the spatial curvature times the 
thermal conductivity k of the medium and the forcing function (S). The thermal conductivity is 
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assumed isotropic within the media, thus, the heat equation model is appropriate for absorbing 
fluid. Heat deposition from an external laser source is the first physical phenomenon to be 
tackled in the treble partial differential system, equations (1-3). 

The force S causing the initial thermal disturbance could be attributed to many sources, 
but we concentrate upon disturbances caused by Q-switched laser pulses. Inspection of the 
forced heat equation shows there are two contributions for rapid, concentrated heat deposition; 
specifically, either the spatial curvature of the deposition is large, meaning the spatial gradient is 
large, or the time derivative of heat deposition is large. In general, the amplitude of the resulting 
pressure wave is directly proportional to the power density of the impinging laser pulse. 

The impinging light is focused on the absorbing media and the absorbed energy raises the 
temperature of the media. Photons entering the media are absorbed according to Beer-Lambert- 
Bouguer's law of absorption, where absorption of energy can be treated omnidirectional; thus, 
the forcing function may be treated as spherically symmetric. The original spatial domain will 
be mapped to the Fourier codomain by employing three-dimensional Fourier transforms (spher- 
ical transforms), where the spherical radial variable (r) is mapped to the codomain variable 
(s), ergo r n- s. Transforms will be represented by either Bracewell's symbolism for a Fourier 
transformation (c) or by the following symbol: J?r3 J6). If the immediate assumption of spher- 
ical symmetry be perceived a limitation, it is possible to break the spatial symmetry at a latter 
point in time by employing the convolution theorem. 

The forcing function S describes the manner in which an impinging laser would interact with 
some ponderable media. It is common to consider the forcing function separable in space and 
time; this greatly facilitates prehending a solution for the partial differential equation. Addition- 
ally, stochastic independence is often assumed for an incident laser beam, given the statistical 
nature of photons. The source function S is separated into the product of two functions, where 
f(t) defines the time dependent behavior of the source, P(r) describes the spatial distribution 
of the laser beam. 

S(r, t) => P /(*)P(r) C P /(*)P(«); r ^ s (4) 

The above equation shows the decomposition of the source function into the spatial and 
temporal components of the incident laser pulse. The three-dimensional Fourier transform of 
the spatial distribution of the incident laser intensity P(s) is multiplied by a suitable constant 
Po to scale the amplitude of the incident laser energy. 

Let the spatial distribution for the light source remain arbitrary for the moment and con- 
centrate upon the particular solution for the forced heat equation. The Fourier transform of the 
forced heat equation, in combination with the source decomposition specified, provides a time 
dependent ordinary differential equation, which can be readily solved Q. The complementary 
solution, commonly referred to as Green's heat kernel, describes the decay of an initial heat dis- 
tribution and is not shown in the solution, because the initial heat profile is assumed regular and 
quiescent. The particular solution describes the contribution of heat attributable to the external 
source acting upon the media over some duration of time, viz.: 

T(M) = QP °J (S) e- a ^ 2t J e a ^ 2t 'f(t') df (5) 

The particular solution of the forced heat equation in the Fourier codomain involves the 
superpositioning of both fundamental solutions to the ordinary differential equation. Integrat- 
ing Green's unbounded heat kernel multiplied by the specific time dependence for the external 
source of light describes the accumulation of heat within the media. The integral is then mul- 
tiplied by Green's heat kernel, the diffusivity a, the magnitude of the laser intensity Po and 
the Fourier transformed spatial distribution of the laser light P(s), and finally, divided by the 
thermal conductivity k. Diffusivity describes the rate at which heat will flow in a media and 
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is defined as the division of the thermal conductivity by both the material density and heat ca- 
pacity. The specific spatial distribution of the laser intensity P(s) will be ignored presently, 
because this term is time independent and therefore may be treated as a constant. 

Regarding the time function f(t'), this describes how the light source is modulated in time, 
but the specific modulation is immaterial. Duhamel's principle can be applied later and con- 
volve the solution with any specified modulation; as a consequence, let the time function f(t') 
be constant over the interval [0, t], which will set the limits on the time integration. Completing 
the time integration leads to a compact solution, viz.'. 



The temperature solution is rather general in nature, the specific spatial distribution of the 
source is multiplied over the kernel solution in Fourier space. The second term (within brackets) 
represents the upper limit of the time integral and is left open-ended for arbitrary time t. In the 
limit of infinite time, the second term vanishes, thus, showing the purpose this term has for 
regulating the time deposited heat within the material. Division by the thermal conductivity n 
scales the magnitude of the resulting temperature. 

Considering the spatial distribution function P(s) a constant, then the solution presented in 
equation (O is essentially a Green's function for the corresponding forced ordinary differential 
equation. To see the influence an incident laser possessing beam characteristics and material 
interaction has on heat deposition, then both the beam width (u>q) and the penetration depth (7) 
must be brought into the solution. The penetration depth is inversely proportional to Lambert's 
absorption for a material and is measured in units of length; equivalently, the penetration depth 
is inversely proportional to the absorption coefficient, i.e. a = 1/7. The beam width represents 
the axial distance for the Full Width Half Maximum (FWHM) of a Gaussian beam. Appropriate 
division by a rational polynomial of the incident Gaussian beam is taken to ensure the Fourier 
transform has unit area, viz.: 



It can be seen that the incident light has a constant power profile, it is most efficient to 
control the intensity by an arbitrary constant rather than parameters used to describe beam 
characteristics. Using equation © to describe the impinging laser light, a specific solution for 
the temperature can be derived, which includes the beam characteristics, viz.'. 



The overall temperature is controlled by the power of the incident light source; Po, mea- 
sured in units of power (Watts). The temperature in the media is mitigated by several factors, 
namely, the geometric dependence (1/r), the beam width (u!q) and the penetration depth (7). 
The algebraic sum of the two error functions represents the cumulative energy absorbed by the 
system and that power is scaled by the width of the incident laser beam, the absorption of energy 
governed by the penetration depth, and the period of time the laser is active. The temperature 
reaches it's maximum as time progresses; this is the particular solution for deposited energy 
from an impinging light source. If diffusion should also be modeled, then the accumulated heat, 
equation dHJ, must be the leading coefficient for Green's heat kernel, which describes the dis- 
sipation of energy within the system. Due to the maximum principle, one cannot represent the 
rise and dissipation of heat on the same timeline; rather, one is forced to segment the two time 
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regimes (8). For energy deposition over large time periods, diffusion becomes the predominant 
process; but, for short time regimes, acoustic wave generation predominates. 

The spatial distribution of the light source determines the spatial relationship temperature 
has with the media; as the width of the incident beam is increased, more of the incident energy 
is distributed over a wider volume and this effectively reduces the magnitude of the temperature 
reached over any period of time. Similarly, the magnitude of the temperature is inverse to the 
penetration depth; hence, the less absorptive the media is to a particular frequency of light, the 
lower the resulting temperature within the media. The contribution of either the beam width 
or the penetration depth can be ignored by simply considering an infinitely absorbing media 
(7 -C 1) and an infinitely tight beam {uiq <C 0); then the deposition of heat is governed solely 
by the product of the diffusivity and the pulse width. For sufficiently short enough pulses, 
measured relative to diffusivity, the ability to realize thermal confinement becomes possible. 
Thermal confinement refers to the process of dumping significant enough energy locally within 
the media, rapidly enough as to cause nonlinear effects. 

Now that the temperature is known, we place attention on the scalar potential. The scalar 
potential function, equation (2), is essentially the forced wave equation and yields a second 
order differential equation after transforming the spatial coordinates to the Fourier codomain; 
once again, assuming spherical symmetry. The forcing function for the scalar potential is equal 
to the time derivative of the temperature function scaled by the thermal expansion coefficient 
(3. The derivative of the forcing function does not involve the spatial variables regardless of 
what domain is considered, either the original domain or the Fourier codomain; thus, the time 
derivative of equation © yields equation (0. 

F(t) ee -p±T(s, t) = -pS^l e -^s)H (9) 
dt k 

The solution for the scalar potential is a sum of two complementary solutions and the par- 
ticular solution; we are only interested in the particular solution. The complementary solution 
provides the behavior of the space itself, also, any constant waves existent within the media prior 
to the action of the forcing function are not of any interest. The particular solution Y(t) involves 
both complementary solutions, y\{t) and U2(t), the forcing function F(t) and the Wronskian 
W(t), viz.: 

With complementary solutions: 

yi (t) = e- 2 ™ vt ,y 2 {t) = e 2 ™ vt (11) 

The Wronskian W(i) describes linear independence for the solution and, by inspection of 
equation ( fT2l . the solution is independent for all time except for the origin of the transform 
variable or a zero velocity wave, i.e. {s, v} = 0. An acoustic wave with a null velocity must be 
rejected. An association with the origin is implicit, because the initial condition is considered 
quiescent; hence, there can be no acoustic wave if the impulse causes a constant disturbance 
throughout space. Because the field is considered constant prior to the action of the forcing 
function, that is, the incident laser light, we must disregard the solution dealing with time zero. 



W{t) = dot 



2/1 (*) y a (t) 
y[(t) y' 2 (t) 



= Anivs = (12) 



The particular solution is generated by considering both the incoming and outgoing forms of 
the acoustic wave. Since the outgoing wave, emanating from a source or disturbance is sought, 
the particular solution Y(t) takes on the form as shown in equation ( fl"3l ). Integrating the forcing 
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function with the incoming wave form describes how the disturbance accumulates over time. 
The time integral is multiplied by the outgoing wave form, describing the transmission of any 
accumulated energy out from the source. 

Now the limits of integration are over the interval of time the light source should cause heat 
deposition within the media, which has already been stated to range from zero to arbitrary time 
t, i.e. [0, t]. Because the quiescent period prior to external perturbation by the incident light 
is inconsequential, as evidenced by the Wronskian, then all that is required is an end point in 
time when the laser pulse is shut off. The restraint for the limits of integration can be succinctly 
represented by implanting a Dirac delta function shifted by the pulse width (to). This is an 
approximation applied to the time integration for the formation of the disturbance and will not 
accurately describe the formation of the acoustic wave itself. The approximation will become 
better as the pulse width in time is shortened. Nevertheless, the chosen method of describing 
the time dependence is applicable no matter the time period for the pulse as long as the range 
for the time variable is truncated accordingly. 

p — 2iTisvt poo 

Y(t) = - / e 2™vt> e - a (2irs) t' § n> _ y dt / (13) 

Airivs J_ 00 

The leading coefficients are not shown in the particular solution Y(t) for brevity, such as 
Po, a, &c. Taking the limit of infinite time, Green's heat kernel diminishes and the thermal con- 
tribution to the resulting wave approaches zero, canceling any acoustic wave; thus, the impulse 
must be chopped or modulated in time to cause an appreciable acoustic wave. This justifies the 
approximation for the integral by truncating to the upper time limit for the laser pulse width. In 
the limit of the pulse width approaching zero, the integral approaches unity, thus, the resulting 
acoustic wave is that from the most fundamental impulse possible, i.e. a Dirac delta function. 
Shifting the placement of this impulse in time approximates the accumulation of deposited heat 
over pulse width. After completing the integration, the solution is a product of a sinusoidal 
wave and an exponential, also, division by the Wronskian, viz.: 

Ap—iirisvit—to) , 
Y(t) = e -«(2-) 2 to (14) 

4ttvs 

The denominator indicates dispersion of the wave, where the magnitude of the resulting 
wave is scaled by the spatial frequency component, thus, the amplitude of the high frequency 
wave component is more swiftly reduced relative to the slower longer period waves. 

It is now a matter of returning the particular solution for the scalar wave equation to the 
original spatial domain; the inverse Fourier transform involves complex integration techniques 
and both real Ji and imaginary 3 parts are shown in equations (Q3) and ( TToT l. 

J sm(2Trsv(t - t )) r _ a{ 2K S ft a \ = 
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The inverse Fourier transform of the imaginary part (3) of the particular solution 
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Where the imaginary error function is defined as such ■ 

erfi(z) = —ievf(iz) 



(16) 
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Both solutions are implicitly valid for the entire real line for both spatial and temporal coor- 
dinates; but, due to the principle of causality, the appropriate range for the spatial coordinates 
should not exceed the distance allowed by the speed of the acoustic wave multiplied by the time 
traveled, viz.: 

{r\r < v(t - t ), r G [0, oo]} , {t\t e [t , oo]} (18) 

There are two solutions generated from equation (TPfl) . each represents a 90° phase differ- 
ence. The sinusoidal solution implies a disturbance arising out of phase with the perturbing 
force causing the disturbance; contrary, the cosinusoidal solution describes a disturbance con- 
current with the perturbation. It is the cosinusoidal solution we are most interested in. Both 
solutions show the amplitude of the resulting acoustic wave is inversely proportional to the 
speed of the wave in the media (v) and the geometric dispersion (1/r). The product between 
diffusivity and pulse width (at ) determines the degree of thermal confinement. The larger the 
diffusivity, the swifter the thermal conductivity of the media, the smaller the amplitude for the 
resulting acoustic wave. Diffusivity restrains the pulse width of the impinging light to enforce 
thermal confinement, where for short enough pulses ablation can occur. In contrast, for pulse 
widths many orders of magnitude larger than the diffusivity, the material allows deposited heat 
to be diffused away, thus, quenching any anticipated non-linear effects. 

Considering the specific parameters attributed to the light source, equation (0, it is a simple 
matter to fold this description into the known inverse Fourier transform solution represented by 
the cosinusoidal part of the wave solution, equation ( TTol i; moreover, this combined solution will 
be signified by the symbol K(r, t) for later reference. 
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Analyzing the solution, the amplitude of the resulting acoustic wave is inversely propor- 
tional to the distance (r), the product (ato), the penetration depth (7), and finally, the radial 
width of the laser pulse (u>q). In general, for either a weakly absorbing material, 7 3> 1, or for 
greatly dispersed beams, uio ^> 1, the resulting amplitude of the acoustic wave is vanishingly 
small; thus, both of these parameters control the ability to cause an acoustic wave of any signif- 
icant amplitude. But, taking the limit of both of these parameters to the other extreme, infinitely 
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absorbing and infinitely confined beam ({7, luq} — > 0), the ability to cause thermal confinement 
is the limiting factor responsible for the amplitude of the acoustic wave. 

The real (5ft) part of the functional, K(r, t), is what is of interest. Keeping in mind the 
constraints causality placed upon the range for all relevant parameters, the functional provides 
a compact, closed-form analytic solution for the scalar potential wave arising from Q-switched 
laser pulses. The solution for the scalar potential ($) can now be constructed, where earlier 
ignored leading coefficients can be reintroduced, viz.'. 

K 

Based upon the known functional for the scalar scalar potential: the velocity v(t), pressure 
pit) and density p(t) of the resulting acoustic wave are all proportional to either the spatial 
gradient (V) or the time derivative (d/ dt) of the scalar potential, viz.: 

velocity : v(t) = -V$ (21) 
d$ 

pressure : pit) = —po—r- (22) 
at 

Pa d$ 

density: pi =-^-7T (23) 
V s at 

3 CLOSING 

A need for a compact, closed-form analytic solution describing laser induced optoacoustic wave 
generation was identified. The most fundamental and efficient solution was sought to serve as 
a kernel for further complications in modeling & simulation and/or theoretical investigations. 
The system of partial differential equations used to describe the absorption of laser energy, sub- 
sequent heat generation, and finally, formation and propagation of an acoustic pulse requires 
multiple integrations; furthermore, it is most efficient to incorporate most of the information 
contained in the treble differential system in as succinct a formulation as possible. A suitable 
solution for the scalar potential enables calculating the velocity, pressure and density of the re- 
sulting acoustic wave. A kernel solution, as the nomenclature implies, serves as a seed or kernel 
for further complication, where specific temporal or spatial constraints can be entertained. 

The analytic solution presented assumes both spherical symmetry in physical space and 
rectangular time dependence. Duhamel's principle allows further refinement of the time de- 
pendence by appropriate convolution of the kernel solution with some suitable time function; 
thereby, enabling temporal modulation of the impinging laser pulse. For the kernel solution 
provided, the beam characteristics attributed to the impinging laser pulse was assumed to be 
spherically symmetric; but, if further complication be desired, convolving the kernel over some 
specific spatial distribution can be performed, where appropriate consideration of the spatial 
basis set is taken. In either case, application of Duhamel's principle or spatial convolution, both 
of these operators should be applied directly to the scalar potential ($) solution and applied 
before any other derivative is applied, specifically, before applying those derivatives indicated 
in equations describing various physical observables, i.e. velocity, pressure. 

The mission of NAMRU-S A/Directed Energy Biomedical Research Department is to inves- 
tigate nonionizing radiation bioeffects and understand the potential risks to human health and 
safety. Q-switched laser pulses can cause damage to the retina, including hemorraghing caused 
by a laser induced optoacoustic pulse. Pulsed laser energy is unique in that the radiant energy 
is very low, but the radiant power of the pulse is very large; thus, the temporal nature of short 
pulses can initiate an array of unexpected physical responses, e.g. optoacoustics, ablation, &c. 
In addition to bioeffects, optoacoustic spectroscopy is a very powerful detection method and 
considerable interest exist to develop a novel method for measuring Specific Absorption Rate 
(SAR) spanning the radio frequency band, specifically, ELF to terahertz. 
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